Asymptotics of eigenvalues for an energy operator of the 

one model of quantum physics. 

E..Yanovich 

St-Petersburg State Polytechnical University , 

Department of Higher Mathematics, 

29 Polytechnicheskaya st., St. Petersburg, 195251 Russia 

E-maih teduard@land.ru 



O 
(N 

>^ 






In this paper we consider eigenvalues asymptotics of the energy operator in the one of the most 
interesting models of quantum physics, describing an interaction between two-level system and 
harmonic oscillator. The energy operator of this model can be reduced to some class of infinite 
Jacobi matrices. Discrete spectrum of this class of operators represents the perturbed spectrum 
of harmonic oscillator. The perturbation is an unbounded operator compact with respect to 
unperturbed one. We use slightly modified Janas-Naboko successive diagonalization approach and 
some new compactness criteria for infinite matrices. Two first terms of eigenvalues asymptotics 
and the estimation of remainder are found. 
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1 Introduction and main results. 

We consider the energy operator of the following form 



H 



huJn 



a^ + hw a^ Si + h\ ((j+ + o-_)( a + a"* 



where a^ , 0"+ , o"„ are the 2x2 matrices of form 



o-z 



1 
-1 



CTa 



1 
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a and a+ are the creation and annihilation operators for harmonic oscillator, A is the interaction 
constant, u is the oscillator frequency, uq is the transition frequency in the two- level system. 
These matrices and operators are satisfies by the following commutative relations 



cr4-,a_ 



a. 



F^,o-^ 



2o-+, [o-^,(J-] 



-2(T_ 



a, a 



1 



In the work pQ it was shown that the hamiltonian of this model is represented by the two 
Jacobi matrices. These matrices have the following general form 



A 



/ ci gVl 

g^A 1 + C2 gV2 

gV2 2 + ci 
gVs 


V 







gV^ 

3 + C2 gVi 

gy/A 4 + ci 



\ 



;i-i) 



/ 



where g, Ci, C2 are real parameters. As it is well known [21 |3] that the matrix A defines a 
selfadjoint operator with simple spectrum and the domain D{A) is dense in the space /2(I^)- 



As the operator A can be considered as relatively compact perturbation of the main diagonal, 
its spectrum is discrete. 

The main goal of this paper is the investigation of the eigenvalues Xn{A) behavior for large 
values of n and fixed other parameters. The result is given by the following asymptotic formula 
(Theorem (13AD ) 



Xn{A) =n-g^ + 



Cl +C2 



O 



1 



1/16 



n 



n — 7- oo {g y^ 0) 



2 Selection of the main component in the asymptotics. 



Let's present the operator A in (11. ip in the form 



A = A, + K^I+'-^^R. 



2 2 

where I is the identical matrix, Aq and R are defined in the following way 

/ gVT ...\ 

g^fl 1 g^f2 

(7^2 2 (7^3 

g^ 3 



(2.1) 



A. 
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R 



/ 1 
0-100 
10 
0-1 
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V 
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The matrix Aq represents so called "shifted oscillator" operator: a^a + g {a + a'^), where a+ 
a are the creation and annihilation operators. If we use its matrix representation, we obtain 
exactly the matrix Aq. 

Eigenvalues problem for the operator Aq has an exact solution. This solution can be obtained 
by the different ways. For example, with the help of Bogolubov's transformation [4j or by using 
continued fractions [3]. This solution has the form 



^oen = /ine„, Hn = n-g, n = 0,l,2,... 



(2.2) 



Here e„ are the normalized eigenvectors of the operator Aq. Its expansion through the basis 
vectors a„ of the matrix representation (II. ip has the form 



Ua„ 



/ J '-^n,'n 



n=0 




f/„,„ = exp{-g''/2} \l— (7™-"L(r-"H^') = ^f-"H^') , (2-3) 



where the matrix elements Un,m of the orthogonal transformation U are defined by Feynman- 
Schwinger's formula [SIE] (see also [S]); Ln{x) are the generalized Chebyshev-Laguerre's poly- 
nomials [TOl 



L^l 



{n + s]\ 



X] 



n\ 



i=0 



X 



ii + s) 



Ct 



n\ 



i\ (n — i)\ ' 



[s>0) 



Lt^H 



-X 



^^i^LilM 



W. 



M) 



and Un (x) are the normalized Laguerre's functions 



coi'Hx) 



n\ 



[n + s]\ 



-x/2^s/2^(.)(^) 



Uj\^ {x) Uj':^\x) dx = 5n,' 



Let's note here that the comleteness of the "shifted oscillator" eigenfunctions (in coordinate 
representation) for complex values of the parameter g is considered in [9] . 

Let's us find the matrix of the operator R in the basis of the operator Aq eigenvectors. 
Denoting the elements of the transformed matrix as Rk,m and taking into account that Rn,m = 

( — lYdr^rn, WC obtalu 



■'n,rm 



Rk,m = {Rem. eu) = {U^RU)k,m = J2^-ir Un,k U^,m (2.4) 



ra=0 



Let's present the matrix elements Un,m as contour integral 

t/„,™ = exp{-//2}.H^"-™-i, /x'"-ifi-l) expj^jrfx, (2.5) 

V nl zm J \x J I a; J 

c 

where C is a unit circle centered in the origin of the complex plane x. 
Substituting (12. 5p in (12. 4p and summing up over n, we find 

Rk,m = exp{-2g^} Vkh^. g^^'^ -1- j j{xr-^ {x'f~^ exp L'' + 1 - -LU rf^; dx' 



c c 



Contour integrals in this expression can be calculated consistently with the help of residues. 
As a result, we obtain 



Rk,m = {-If exp{-2/} J— {2gr-' Y.CI{-1 



k ' ' -,2\i 



"^' fn^\m-k X^r^if ,M (V 



i=0 

(if the expression {i + m — k) in the denominator of the last sum is negative, the corresponding 
term must be considered equal to zero). 

Comparing this expression with (12. 3p . we obtain 

R,,.^ = {-If Jr'\^9') (2.6) 

In spite of seeming asymmetry, the matrix Rk,m is symmetric {Rk,m = Rm,k)- It can be 
easily verified with the help of the known property of the generalized Chebyshev-Laguerre's 
polynomials [10] 

Li"\x) = {-xr^^^^Lts{^) 
Using the asymtotics of the generalized Chebyshev-Laguerre's polynomials [10] 

L^(x) = TT-^/" n^/2-1/4 3,-./2-i/4 gX/2 |cos(2v/^ - s7r/2 - 7r/4) + 0{n~^'^)] , n ^ oo , (2.7) 

we find 

lim ^„,„+p = , Vp G Z (2.8) 

n— ^oo 

In what follows we will need the following result [TT] 

Lemma 2.1 (J.Janas-S.Naboko) Let D he a selfadjoint operator in a Hilbert space H with 
simple discrete spectrum (Dcn = f^'nen), where {e„} is an orthonormal basis of eigenvectors in 
H and /i„ are simple eigenvalues (/i„ — )■ ooj, ordered by |/ij| < |/ij+i|. Assume that |/ij — /i^l > 
Co > 0, Vi 7^ A;. If R is a compact operator in H then the eigenvalues \n{T) of the operator 
T = D + R (with discrete spectrum too) become simple for large values of n and satisfy to the 
asymptotic formula 

A„(T)=^„ + 0(||i?*e„||), n^oo, (2.9) 

where R* is the adjoint operator with respect to R. 



The matrix R and the matrix similar to it i? = U^RU, represents bounded noncompact 
operator (projector) since R^ = I. Therefore we can't apply here at once the Lemma fl2.ip . 
Let's prove the following theorem: 

Theorem 2.1 Let D be a self adjoint operator in a Hilbert space H with eigenvalues /i„ = n, 
(n = 0,1,2, .. .) and complete system of corresponding eigenvectors in H . Let R be a bounded, 
selfadjoint, noncompact operator and its matrix Rn,k in the basis of operator D eigenvectors 
satisfy to the condition 

lim Rn,n+p = , Vp G Z (2.10) 

Then the eigenvalues Xn{T) of the operator T = D + R (having a discrete spectrum too) become 
simple for large values of n and satisfy to the following asymptotic estimation 

Xn{T) = n + Rn,n + O (sn) , n ^ oo , (2.11) 

where 



I D 12 

and Sn ^ at n -^ oo. 

For the proof of this theorem we need the following compactness criteria for infinite matrices. 

Lemma 2.2 Let V be a bounded, noncompact operator in a Hilbert space H . Let its matrix 
Vij (i,j = 0,1, . . .) in some orthonormal basis satisfy to the condition 



limV;,„+p = 0, VpeZ (2.12) 

n— >oo 

Let b = {&i}^_oo is an arbitrary l2-sequence 

i=— oo 

Then the operator K with matrix Kij = bi^jVij (i,j = 0,1, . . .) is compact in H . 

Proof. Let's show at first that the operator K is bounded. For that we need to prove the 
estimates [12] 



j=0 i=0 

where A is a constant independent of i and j {\\K\\ < A). Using Cauchy's inequality, we have 



oo 



1/2 / oo \ 1/2 



j=0 j=0 \j=0 J \j=0 J 

Due to (I2.13P the first estimate in (I2.14p is fulfilled. By the same way the validity of the 
second estimatein (I2.14p is established. Thus the operator K is compact. Let's prove now its 
compactness. 

Let's define the cut-off function fe*^"^ = {&j }i^_oo of the sequence {&«} 

,(n) r 0, \i\ > n 
bi, \%\<n 



Let's define the sequence of operators K^"''^ by tlie formula KI"- = hiJ-Vij. It follows from 
this definition and from fl2.12p that K'^'^'> is a compact operator for arbitrary n. We have further 
so as at the proof of (12.141) 

||i^_i^H|| < ||6-6(")|| .||V|| 

Terefore 

and K is compact as a limit by norm of compact operators [T^ [T3] . 

Proof of the theorem ( 12. Ih . Let's associate to each operator the matrix in the basis of 
the operator D eigevectors, remaining the same notations. 

Following the main ideas of the work ^lj let's show that there exist such anti-hermitian 
operator K {K* = —K) that 

{I + K)T-Di{I + K) = B, (2.15) 

where B is compact operator and Di = D + diag{i?„^„}. (So Di is the diagonal matrix with 
elements (-Di)„,„ = n + Rn^n)- Suppose that such operator K have found. The condition (I2.15P 
means that 

T = (/ + K)-\D, + B{I + K)-^){I + K) 

(The existence of the inverse operator {I + K)~^ follows from the anti-hermitianess of K). That 
is the operators T and Di + B{I + K)~^ are similar and have the same spectrum. But the 
operator B is compact and the eigenvalues of Di due to (I2.10p satisfy to the requirements of 
Lemma (12. ip . Applying the Lemma (12.10 we obtain 

\n{T) = n + Rn,n + 0{\\B*en\\) (2.16) 

Therefore for the proof of the theorem we should establish the existence of such the operator K 
and find the matrix of the compact operator B. Substituting the expressions for the matrices 
T Di in (I2.15P we obtain 

(/ + K)T -Di{I + K) = Ri~ [D, K]+KR- diag{i?„,„} K , (2.17) 

where [■, ■] is the commutator and Ri is the matrix of the operator R with zero main diagonal 

{Ri = R~diag{Rn,n}). 

This expression will be the matrix of compact operator if we can find such an compact 
operator K that the condition [D, K] = Ri is valid, or in matrix form: Kij {i — j) = (i?i)j j. It 
follows from that 

Ki,j = -^, t^j; Ki,i = 0, 1 = 0,1,- ■■ (2.18) 

As the operator R is selfadjpint the corresponding to the matrix (12.180 operator K is anti- 
hermitian. Its compactness follows from the Lemma (12. 2p . Actually, if we choose the sequence 
{bi} as {1/^} {i 7^ 0) then from (12.100 it follows that all conditions of the Lemma (12. 2 p are 
fulfilled. 

Now from (12.170 we find the form of the compact operator B: 

B = KR-diag{Rn,n}K, 
and taking into account (12.180 we obtain 



0(||B-e„||)=0(||A-e„||) = 



:\ 



I f? |2 

(n - A;)2 



E 



Substituting this estimate to (12.161) . we obtain the formula fl2.1ip . The theorem is proved. 



Due to (12.81) . the condition fl2.10p of the Theorem (12.11) is fulfilled. Hence, applying the 
Theorem (12. ip and taking into account ( 12. 2p and ( 12. ip we have the following result 



Xn{A) =n-g^ + 



Cl +C2 



+ O (Sn) , n — )■ oo , 



where 



IP 12 
l-'^fcjnl 



^S("-'=)^ \ 






in - kY 



(2.19) 



3 Estimation of the remainder. 

To estimate the decreasing rate of the sequence s„ we should have another estimation for the 
Laguerre's functions ujn{x) than the estimation following from (12. 7p (in (12. 7p the parameter s 
is fixed). We could not find this result among known one and therefore we give here not only 
the formulation but the proof of it. 



Lemma 3.1 Suppose that x > 0, s G Z^. Then the following estimate for the Bessel functions 
Js{x) is valid 



\Js{x)\<2J— (l + - 
V vrx V X 



Proof. Let's use known representation 




Js{x) = \l — {P{x, s) cos(x — stt/2 — 7r/4) — Q{x, s) sin(x — sn/2 — 7r/4)) 



where 



^ ' 2r(s + l/2) J U 2a; 





s-l/2 



1 - 



2i 



s-l/2' 



du 



Q{x,s) 



2ir(s + l/2) 



e-"M^-i/2 



. X s-1/2 

2x) 



■ N S-1/2' 



It is evident that 




Js{x)\<\—{\P{x,s)\ + \Q{x,s)\) 



TTX 



(3.1) 



and everything reduces to the estimation of the integrals P{x, s) and Q{x, s). Let us consider 
the integral for P{x, s). The estimation for Q{x, s) is the same. At s = we have \P{x, s)| < 1, 
|Q(x, s)| < 1 and the estimation (13. ip gives the required inequality. Suppose that s & N. In 
this case we have 



\P(x,s)\< 



ns 



+ 1/2) J \ 2x) 



u \«-i/2 



du 



Expanding the binomial in this integral in the series on u/2x 

fc=i 

^ is-l/2).....{s-l/2-{p-l)) ^ ,, ^uy^ ^ ^ 

k\ \2x/ 

and putting p = s we have (1 + oy^P^^/'^ < 1 and therefore 

fc=l 

Integrating by terms we obtain 

IP^, ,^|^1 I y^ r(. + fc + l/2) (.-l/2)-...-(.-l/2-(fc-l)) 1 

in^,sj|<i + ^ r(s + i/2) fc! (2x)^'^ 

For Q(x, s) the same estimate is valid and the formula (13 .ip leads again to the required inequal- 
ity. Lemma is proved. 



Lemma 3.2 If x > 0; n,s E Z^ and s^^ < n then 

( 

^+ 1)1/4 ' 



^i^nx)|<r^^, (3.2) 



where the constant C{x) depends on x only. 

Proof. Let's use Laguerre's functions integral representation trough the Bessel func- 
tions [lo] 



e^/2 



= I e"* t"+t Js{2Vtx) dt, n,seZ^ 



^Jn\ (ri + s)! 
Let us split this integral into two one 

px/2 / /• 

uj^:\x) = , \ + 

^ ' Jn\{n^s)\ \J 

^ \0 to 

where to > is an arbitrary now. 

For estimation of the first integral let's use the known inequality [H] 

\Js{x)\ < 1 , s E Z^ , xER, 

For estimation of the second integral we use more precise estimate from Lemma (13.1 p 

Jsix)<2\ — 1 + - <2e\ — , x>s^ 
V vrx V x/ \ nx 



Putting to = s'^ /Ax (so that at t > to one can use the last estimate) we have 

" to 



^):k^)\ < 



,x/2 



< 



^Jn\ {n + s)\ 

gx/2 



e-'e+2dt + 



2e 



/ 



TTVX 



e-*t"+f-i/^rft 



to 



< 



n\[n + s 

x/2 



to max<ie"*t"+5 K + 



2e 



/ 



T{n + s/2 + 3/4) 



vrvx 



^4 g-(n+2)(^ + 5/2)"+2 2e r(n + s/2 + 3/4) 

4x 



^Jn\{n + s)\ A/TTv/i A/ri! (ra + s)! 

Using known inequahties for Gamma-function following from Stirling's formula 

Ci z'^^l^ e~' < T{z) < C2 z'-^'^ e~^ z > <5 > , 

where Ci, C2 are some constants independent of z, let's estimate each term in square brackets. 
We have 



2e r(n + s/2 + 3/4) 



a/^a/x A/n!(?2 + s)! 



<C7(x) 



[n 



+ s/2 + 3/4)"+^/2+V4 



+ l)n/2+l/4 (^ + s + 2)^/2+^2+1/4 



< 



[n 



^ C{x) (1+2(7^)" ^ C{x) 



1 + 



n/2 



^n+ 1)1/4 (l + _^)n/2 - (^+1)1/4 ^- ' 4(n+l)2 

At last, if we put n> s^ then 

2e r(n + s/2 + 3/4) C(x) 



Cix 
< ; , , . exp 



(n + 1)1/4 - ^ 1^8(^ + 1) 



\/tQx 
Similarly, we can estimate the second term 



^n\{n + s)\ ~ (n + 1)1/4' 



[n> s' 



(3.3) 



4 p-("+f) 



s e 



f)(n + s/2)"+S 



< 



C{x 



4a; ^n\ {n + s)\ 



<C{x) 



s^ {n + s/2)"+^/2 



'n + l)"/2+l/4 (n + S + l)n/2+s/2+l/4 



< 



:i + ICTTTT) )" ^ C7(x) 



(n + 1)1/4 ^n + s + 1)1/4 (2 + ;^)"/2 - (n + 1)1/^ (n + s + l)i/^ 

If ra > 5!^ then 

a6\ 



exp 



n + 1 



4x y/n! (n + s)! 

From (Q, (I33D it follows that 



s4 g-(n+|) ^^ ^ s/2)"+l C(a;) 

^ ' +1)1/4 ' (^ ^ ^ 



[n 



(3.4) 



K-'wiSj;^. ( 



n > s 



16n 



q.e.d. 



Having the estimate (13. 2p and the orthogonality condition of the transformation U: 



5^if/fe,nr - Yl 



\^k 



fc=0 



fc=0 



Z^\^n-k\ + / J^« I 



1 , Vn G Z4 



(3.5) 



fc=0 



k=l 



\u 



(n-fc), _ I (fc-n) 



one can give the estimate of the remainder which is defined by the sum 



(n-fc)|2 n I (k) ,2 00 , (fc)|2 



E l^fc I ^ V^ \^n-k\ , V^ 

(r, - i-^2 Z^ 1,2 "^ Z^ 



W^i 



fc^n ^ ^ fc=l 



fe=l 



A;2 



(3.6) 



From Lemma fl3.2l) it follows that 



\^n—k\ — 



C 



(n- A; + 1)1/2 



n-k>k'^ (n > A;^^ + A;) 



a;(^)p< 



C 



+ 1)1/2 



m 



n>k^'> (A;<ni/i^) 



1/16N 



Let A;„ > be a maximal nonnegative integer of k, satisfying to the equation n > k^^ + k. 
It is evident that k„ < n^/^^. Hence 



(k) 1 2 ^ 

l^n-fcl ^ (n- 721/8 + 1)1/2 

' " ' - (n+ 1)1/2 



rC "^ fvri 



fv *^ rC*n 



Let's present the sum (13. 6 p in the form 



(ra-fe)|2 

('n. — 

ky^n 



(n - A;)2 



(fc) |2 fcn I (k),2 



E%^ + E 



Wn 



fc=l 



fc=l 



A;2 



+ 



" I, ,(fc) 12 oo I (A;)|2 



E 

_K — fi'Tl. I" J- 



W, 



n—k\ 






/C /v)T,H-i 



A;2 



Due to last inequalities we have 



(fc) |2 fc„ I (fc)|2 



E%^ + E 



w« 



.fc=l A:=l 

Since A;„ ~ rai/i^, we have using (13.51) 

"' \, ,W |2 oo I (fc)|2 



< 



(n - ^1/8 + 1)1/2 Z^ p - ^ 1^;^ 



E %^+ E 



UJn 



_/c — /c^iH-i 



K — K-n-^-l 



A;2 



< 



< 



{kn + I] 






o 



(A;„ + l)2 V'^ 



./C Kyi + l 

1 



< 



.1/ 



Combining both estimates, we obtain 



E 

ky^n 



\U 



(ra-fe)|2 



(n - A;)^ 



oi4.Uo^ 1 



n 



1/2 



n^ 



O 



n^ 



Taking into account the formula (12.191) . we come to the following main result 

Theorem 3.1 The eigenvalues A„(y4) of the operator A U.l\) at g ^ have the following 
asymptotics 



XniA) =n-g^ + 



C1 + C2 



+ 



1 



n 



1/16 / ' 



n — > 00 
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